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We consider the implications of the Color Glass Conden- 
sate for the central region of p + A collisions. We compute 
the k± distribution of radiated gluons and their rapidity dis- 
tribution dN/dy analytically, both in the perturbative regime 
and in the region between the two saturation momenta. We 
find an analytic expression for the number of produced gluons 
which is valid when the saturation momentum of the proton 
is much less than that of the nucleus. We discuss the scal- 
ing of the produced multiplicity with A. We show that the 
slope of the rapidity density dN/dy provides an experimental 
measure for the renormalization-group evolution of the color 
charge density of the Color Glass Condensate (CGC) . We also 
argue that these results are easily generalized to collisions of 
nuclei of different A at central rapidity, or with the same A 
but at a rapidity far from the central region. 



I. INTRODUCTION 

The color field of a strongly Lorentz boosted hadron 
can be described as a classical color field Q, so long as 
one has a high enough density of gluons such that the 
field modes have very large occupation numbers. The 
typical transverse momentum scale for which the field 
modes have large occupation number will be called Q s , 
the saturation momentum. Seen with a resolution scale 
of Q s , the collision of two hadrons (say pions, protons, or 
nuclei) at very high energy can be viewed as two (highly 
Lorentz contracted) sources of color charge propagating 
along the light-cone. Renormalization group evolution in 
rapidity leads to a longitudinal extension of the source. 
That is, the charge distribution is spread out on a scale 
given by the characteristic longitudinal momentum of the 
"hard" particles which generate the source of color charge 
entering the Yang-Mills equation for the "soft" modes. 

The field in front of and behind each "sheet" of charge 
is a pure gauge . The color electric and magnetic fields 
associated with these pure gauge vector potentials vanish, 
except in the sheet where the vector potential is discon- 
tinuous (on a scale larger than the longitudinal spread of 
the color charge source). When the two sheets collide, 
corresponding to the tip of the light-cone, the two charge 
sheets interact. This produces radiation in the forward 
light cone. 



The point of our paper is to compute this radiation for 
collisions of particles with different saturation momen- 
tum scales. This problem turns out to be more tractable 
than that of collisions of two particles with equal satu- 
ration scales. For example, to compute the production 
of particles in the central region of equal A nuclear col- 
lisions, one must perform intensive numerical computa- 
tions |Q] . If one collides protons with nuclei at very high 
energies and studies the central region of particle produc- 
tion, there are two scales, the saturation momentum of 
the proton and that of the nucleus. In the limit where 
Aqcd < QP™ ton < Qf, we shall see that the problem 
simplifies, and one can obtain analytic results for quan- 
tities such as the total multiplicity density of gluons at 
zero rapidity. 

The saturation momentum squared is proportional to 
the total number of gluons in the hadron wavefunction at 
rapidities larger than that at which we compute the pro- 
duction of particles. One could introduce an asymmetry 
in the saturation scales by considering equal A nuclear 
collisions far from the central rapidity region. Alterna- 
tively, one could consider collisions of different nuclei, or 
various combinations of the above. In this sense, the pro- 
ton in the p + A scattering case which we consider should 
be thought of as a generic acronym for asymmetric nu- 
clear collisions in either baryon number or rapidity. 

We shall specifically consider the situation where the 
source propagating along the x + = (t + z)/y/2 axis is 
much weaker than that propagating along the x~ = 
(t — z)/y2 axis. In such a case, the saturation momen- 
tum scale Qs^ on which source one can be viewed as a 
classical field is smaller than the corresponding scale for 

(2) 

source two, Q s . This fact has a very interesting con- 
sequence. Namely, we expect three distinct regions in 
transverse momentum. At large transverse momentum, 

(2) 

k± > Qs , both fields are weak. Thus, perturbation the- 
ory should be a valid approximation in this regime [^^) • 
On the other hand, for Q s 2 ^ > k± > Q S X \ the field one 
is weak, and can be treated perturbatively; but field two 
is "saturated", that is, the field strength F^ v = 0(1/ g) 
has attained maximum strength GL^[, and is in the non- 
linear regime. In that regime of transverse momentum, 
field two can not be treated as a small perturbation, even 

(2) (2) 

if > Aqcd and the coupling a s (Q s ) <C 1 is weak. 
Those non-linearities modify the transverse momentum 



f 



distribution of radiation produced due to the interaction. 
Our goal here is to compute the distribution in the in- 
termediate regime Qs > k± > Qi. Finally, at an 
even smaller transverse momentum < Qi X \ both fields 
are strong. In this region we expect a flat k± distribu- 
tion, up to logarithms of k\. However, we can presently 
not compute the distribution in that region analytically, 
but it has been obtained numerically m . 



t 




FIG. 1. The solutions of the Yang-Mills equations in the 
various parts of the light-cone. The charge distributions prop- 
agate along the x~ , x + axes. In the space-like regions behind 
the charge distributions the fields are just gauge transforma- 
tions of vacuum fields, rotated by the respective charge den- 
sities of the sources. In the forward light-cone, the field at 
time — » oo is given by gauge rotated plane wave solutions (3 
and (5\ 

The solution of the nonabelian Yang-Mills equations 
is illustrated in Fig. |l|. The two color-charge distribu- 
tions propagate along the x~ , x + axes. The fields A 1 in 
the space-like regions behind them are just gauge rotated 
vacuum fields, and A = there. 

For an abelian gauge group (electrodynamics) the field 
in the forward light-cone is just the sum of the two 
pure-gauge fields behind the propagating charge distri- 
butions. It is also just a gauge rotated vacuum field, and 
so no radiation occurs (if recoil is neglected). For a non- 
abelian gauge group (chromodynamics) the sum of two 
pure gauges is not a pure gauge, and radiation occurs 
at the classical level, even when recoil is neglected. At 
asymptotic times, the field in the forward light-cone must 
be given by gauge rotated plane waves. In a leading or- 
der perturbative computation |^-|5| those gauge rotations 
can be expanded to first order in the gauge potentials. 
However, in order to reach into the non-linear "saturation 
regime" we must account for the interaction of the radia- 
tion field with the fields of the color-charge distributions 
on the light-cone to all orders. A numerical approach 
to this problem has been used in ref. for collisions of 



equal-size nuclei, and at midrapidity. For current-nucleus 
interactions (Deep Inelastic Scattering off large nuclei) 
the distribution function of produced gluons in the frag- 
mentation region has been obtained analytically (via a 
diagrammatic approach) in 0, where the authors also 
discuss the generalization of their result to p + A colli- 
sions. In the central rapidity region (the region z -C t in 
Fig. |l|) one should also account for the renormalization 
group (RG) evolution of the CGC color charge density 
per unit transverse area [HJ9). The purpose of this pa- 
per is to derive analytically an explicit expression for the 
transverse momentum and rapidity distribution of pro- 
duced gluons in A\ + Ai collisions at high energy, valid 
at all rapidities where the CGC color charge density in- 
cluding RG evolution is much larger for the source A2 
than for A\. Our explicit result for dN/dk\dy shows 
that the transverse momentum distribution is modified 
from a ~ behavior in the perturbative regime (high 
k±) to ~ m ^ ne region where k± is between the 

saturation scales for the two sources. Furthermore, we 
show that the slope of the multiplicity per unit of rapid- 
ity, dN/dy, provides an experimental measure for the RG 
evolution of the CGC color charge density. 

This article is organized as follows. In section || we de- 
rive the transverse momentum and rapidity distribution 
of the radiated gluons to all orders in the field of the large 
nucleus. We do this by solving the Yang-Mills equations 
with the appropriate boundary conditions. This section 
is somewhat technical and can be skipped by readers in- 
terested only in the results relevant for phenomenology 



In section III we discuss the most important features of 



the radiation spectrum in the perturbative regime (high 
transverse momentum) and in the "saturation regime" , 
including the A-scaling and the evolution in rapidity. We 
outline possible ways of measuring experimentally the 
RG evolution of the color charge density of the CGC. We 



summarize in section [V 



II. 



THE DISTRIBUTION OF PRODUCED 
GLUONS 



In this section, we shall first solve the Yang-Mills equa- 
tions in coordinate space. We assume that in the forward 
light-cone the vector potential depends only on the trans- 
verse coordinate x± and on proper time, r = \/t 2 — z 2 = 
V2x + x~ (which is invariant under longitudinal Lorentz 
boosts), but not on rapidity y = log(a: + /a: _ )/2. When 
performing the path integral over the "hard" source for 
the classical color field in eq. (^) we shall explicitly con- 
sider the dependence on rapidity. 

In the space-like regions the transverse fields are pure 
gauges 0, 



= --u„ 



)d i ul{x A 
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= -L e -i9* m {xA.) d i e i9* m {*x) ( m = 1,2) . (1) 



a l 3 (T = 0, x±) = a\(x±) + a 2 (x±) 



The fields $ m satisfy 

- V\<&m = gp m (x±) 



(2) 



We take the distribution of the sources of the gluon color 
field for each nucleus as a Gaussian according to the 
McLerran-Venugopalan model, 



J Vp\T>pi exp(-Fi[pi] - F 2 [p 2 ]) , 



where 



Fi\Pi] = dyd x± trp i (x±,y)/fi i {y) 



(3) 



(4) 



The quantity p, 2 (y) is the color charge squared per unit 
rapidity and per unit transverse area scaled by N£ — 1. 
It can be related to the gluon distribution function with 
known coefficients, as shown in Ref. It will turn out 
that the radiation distribution depends only on integrals 
over /i 2 (y), i.e. the total color charge squared from ra- 
pidities greater than that at which we are interested in 
computation. 

The field $i will be assumed to be weak such that the 
exponentials in eq. (0) can be expanded to leading order, 



(5) 



In the forward light-cone we write the transverse and ± 
components of the gauge field as 



A 1 (t,x ± ) = 4(t,xl) , 
^4 ± (t, x±) = ±x ± a(r, xj_) , 

corresponding to the gauge condition 

x + A~ +x~A + = . 

Thus, our ansatz for the gauge fields is 



(6) 
(7) 



(8) 



A'(t,x ± ) = al{T,x ± )Q( X -)e{x + ) 
+ a\{x ± )G( X -)G(-x + ) 
+ al(x ± )0(-x-)0(x + ) , (9) 
A ± {t,x±) = ±x ± a(T,x 1 _)0{x-)®(x+) . (10) 

Next, we determine the boundary conditions for 
x~,x + — > 0. In that limit, 

[D+,F +l ] + [L>_,F- J ] = 25(x-)6(x+) (a* - a\ - a 2 ) . 

(11) 

(The contribution from [Dj,F^] is not singular at x + — 
x~ = 0.) For this term to vanish identically we must 
satisfy the boundary condition 



as found before in 

Using (|l2|), the equation 

[D^F»+] =J+ = 5(x-)gp 1 (x ± ) , 



(12) 



(13) 



where pi is the charge density per transverse area in nu- 
cleus 1, gives for x + , x~ — > 

8(x~) {2aQ(x + ) + d t a[ -Q(x + )ig [a\,a l 2 ] } 
= S(x-)d i a\. (14) 

That requires the matching condition (3|Q| 

(15) 



y(t = 0,x ± ) = — [a\(x±),a 2 (x±)] 



We now determine the solution in the forward light- 
cone, x+,x~ > 0. [D^F^] = becomes § 



-d T T 3 d T a~ [D\[D u a]] = , (16) 
[D t ,d T al] +i 9 T[a,d T a] = , (17) 



-d T rd T al - igr 2 [a, [D\a\] - [D\F jl ] = . (18) 

We assume that the field of the second nucleus is much 
stronger than the radiation field, and so linearize the 
equations of motion in a. (Note that a —> if source 
one becomes arbitrarily weak, as no radiation occurs in 
the single nucleus case^ That amounts to dropping the 
second terms in eqs. ( |17|Jl^ ). We perform a gauge rota- 
tion 



a(r,x±) = U 2 (x±_)0(t, x±)U£(x±) , 



(19) 



a\{r, x±) = U 2 ( Xl _) I P(t, x±) - -d l j U^x) , (20) 

with U 2 = exp(— ig§ 2 ) as defined in (0). Then [D % ,-] 
becomes the ordinary derivative d l up to corrections of 
order 0((3 l ) which do not show up in the linearized equa- 
tions of motion, 



\d T T Z d T (i - d l d l (3 = 

T 

dr&p 1 = 

-d T Td T f3 l - (d k d k 5 ij - d l d 3 ) ^ = 



(21) 
(22) 
(23) 

Gauge rotating the boundary conditions ( |12| , |15p gives 

I3\t = 0,ar x ) = U^{x ± )a\{x 1 _)U 2 {x ± ) , (24) 
(3{t = 0,x ± ) = 

l ^Ul(x x ) [a\(x ± ),a 2 (x ± )] U 2 (x x ) . (25) 

From (p2[), 9i/3 J (r, x±) is time independent. We can thus 
write 



3 



(3\t, Xi _) = e u d l xir, Xl _) + 5 l A(xJ 



(26) 



where e l1 is the Levi-Civita tensor in two dimensions. 
The first term contributes to the curl while the second 
contributes to the divergence of (3 l . This ansatz for (3 % 
makes ( p2[) an identity. The equations of motion ( pl[]2^ ) 
now read 



-^d T T 3 d T {] - d l d l /3 = 
-d T rd T x - d l di\ = 



(27) 
(28) 



where x = —d J djX- The boundary condition for \ can 
be obtained by noting that x — (3 l , 



X (t = 0,x x ) = e ij d j Ul(xx)a\(x x )U 2 (x x ) 



(29) 



The equations of motion ( [27|]28] ) are solved by a super- 
position of Bessel functions, 



0{t,Xx) = 
x{t,xx) = 



d 2 k 



(2tt)2 

d 2 k 



h{kx) — Ji(wt) , (30) 

LOT 



(27T) 



b 2 (kx)Jo(uT) . 



(31) 



The functions 6i(fcj_), b 2 (Asj.) are determined in coordi- 
nate space by the boundary conditions for the fields at 
r = 0: 



b 2 (xj_) = 
h(xx) = 



C± „»fej 



(27T) 2 
d 2 k 



-h(kx)=x(r = 0,xx) , (32) 



(2tt) 



± e <fc - L -*- L &i(kx) = 2/3(T = 0,»x) • (33) 



This follows from the expansion of J v (x) for small x: 
J„(a;) ~ {x/2) v jv\. Asymptotically, for x — * 00, the 
Bessel functions are J v (x) ~ yj2/irx cos(x — ^tt/2 — 7t/4). 
Also, for time r — > 00 we assume free fields, cl> = |fcj_|. 
Then, comparing the solutions (^,^l|) for r — * 00 to 
plane waves ||, 



d 2 k. 



1 



(3{t — > 00, a;x) = 



£ l (r -> 00, x_l) - , — 

7 (27r)^ V2wr w 

x {a 2 e^-^-^ T + c.c.} , 
yields for purely imaginary b 2 (k±) and real b\{k±): 

oxffcj.) = - 7 l—b 1 (kx)e 3 ^ 4 
y/irkx 

l 9 r 3Tri/A I ^2, -ikj_-x± 



(2tt) 2 V2^3 
x { ai {kx)e lk ^ x ^- lu)T +c.c.) , (34) 
d 2 fc, 1 e u k\ 



(35) 



V^X 

xUl(xx) [a{(xx),a 2 (xx)] U 2 {x ± ) , (36) 



a 2 (k±_) = 



1 



i 



-b 2 {kx)e' 



r/4 



Xe tf ^C^(arx)ai(i±)Cr 2 (a5x) . 



(37) 



To simplify a\, recall from (Q) that a\ = U 2 {—\/ig)d l U\. 
Therefore, 

igU\ [a\,a 2 ] U 2 = d l U%a\U 2 ~ U\ (d l a\) U 2 

= a{ l d l ult a U 2 . (38) 

Let us evaluate tr|a2| 2 first. Squaring the amplitude and 
taking the trace yields 



tr « 2 



1 



Trk 2 



d 2 x ± d 2 



±a z±e 



-ik±-(xj_-z±) 



x e^e kl did l z tr(A* (x ± , y) A k (z ± , y)) $1 , 



(39) 



Here, y denotes the rapidity. The averaging is with re- 
spect to the gauge potentials $1 and $2, assuming a 
Gaussian weight QS| : 



(0)$ = Jv<PO(<S>) 



x exp 



dy / d a: x 2 2( ^p- 



(40) 



When averaging over $2, the j/'-integral extends from 
— 00 (or some large negative rapidity beyond which the 
source vanishes) to the rapidity of the produced gluons, 
y. Vice versa, when averaging over $1 it goes from y to 
+00. 

We now have to compute the correlation function 

(41) 

with 

A i (x ± ,y) = Ul(x ± ,y)(d i $ 1 (xx,y))U 2 (x ± ,y) . (42) 



The average over $1 in eqs. (|3£ J41|) can be performed 
right away. From (^) we have |)| 

did k z {^l{x u y)^\{z u y))^ 



g 2 5 ab / dy>f(y')« 7 (ux) 
Jy 

g 2 5 ah Xl {y)dld k 1 {u 1 _) , 



(43) 



with u± = x± — z±. Also, we defined the total charge 
squared at rapidity y induced by the source from rapidi- 
ties [y, 00] (not to be confused with the auxilliary fields 
X, X used above in intermediate steps of the calculation), 



Xi(y) = / dy' n\(y') 
Jy 



(44) 
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The tadpole-subtracted propagator is 
1 



7(^-0 



— x\ loga^A 2 ^ 



(45) 



Also, in ((43|) we assumed slow variation of [i\ over the 
relevant transverse scales, and so neglect derivatives of 
it. 

We are left with 

tr(ut(x ± ,y)t a U 2 (x ± ,y)Ul(z ± ,y)t a U 2 (z ± ,y))t, 2 . (46) 

The most efficient way to evaluate this expression is to 
note that 



U}{x ± ,y)1»U2{x ±t yj) = (V) U^(x ± ) 

/a/3 \ / a/3 



( 

(V) ^ (v exp |" dy' T b $^, y')) ) • (47) 
The path ordered exponential can be expanded as 
1 + ig [ dy'T b , a $ b 2 (x u y') 

J — OO 

J —oo J — oc 

+ ■■■ (48) 

According to ( f46| ) we have to multiply two such expres- 
sions, one at x± and the other at z±. The zeroth or- 
der is of course trivial. The contribution to Q{g 2 ) arises 
from the product of the two O(g) terms in ( |48| ) because 
tadpoles only enter via a subtraction of the propagator, 
j(u±), at u± = §. We find 

( l9 )i\ y '[\y'T b al Xi(^(^,y')^(^y r ) 



= g 2 N c S bb ' dy' dy' ($>(x x , y')$l(z x , y')), 

J — oo J — oo 

= g i N c 5 bb ' 1 (u ± ) Cdy'^y') 



g*N c 6 bb 'j(u ± )x2(y) 



(49) 



Analogously to the definition of Xi(y) above, x^iy) de- 
notes the total charge squared at rapidity y induced by 
the source from rapidities [— oo, y], 



(50) 



X 2(y)= 



Next, we multiply two terms of 0(g 2 ) from eq. (48). 
Again, besides a subtraction at u± — tadpole diagrams 
can be disregarded, and so this is the only contribution 
to that order. 

/y ry' ry ry' 
dy' / dy" / dy' / dy" 
-oo J — oc J — oo J — oo 

v rpb rpC rpb rpC 

x 1 a' d 1 da 1 a' d' 1 d' a 



x (*\{x x , y')mx ± ,y")$« {z x ,tf)*{ {z^y") 



(51) 



We can now contract ^\ with $2 (and accordingly 
with $2 )i or we can contract $2 with $2 (and accord- 
ingly $2 w ith <&| ). However, the latter is zero because 
of the ordering in rapidity. Thus, we obtain 



g*N 2 cl 2 { Ul _) f V dy' f y dy"^(y')4(y") 



— 00 J — 00 



W(«±)^(w) 



(52) 



One can repeat the above steps to any order. Resumming 
the series and summing over the one remaining adjoint 
color index we find for eq. (Eq) 



N z 



1 



exp{g 4 N c 7 (u±)x2(y)} 

The correlation function (^d]) reads 
N 2 -l 



(53) 



-g 2 Xi{y)[dld k zl {u±)} 
x exp {g 4 N c j(u ± )x2(y)} ■ 



(54) 



From @, e ij € kl dld l z acts on (g). The direct prod- 
uct with d x d k j(u±) gives zero, such that effectively 
e^e kl d x d l z acts on the exponential only. 

Using (B8|) in (pfl) one derives a very similar result for 



tr|ai| 2 , with the replacement e lj e kl d x d l z — ► S lj S kl dld l z , 
and where again these derivatives act on the exponential 
only. In total we obtain 



tr(|ai 



N 2 



1 



«2 



2nki 



d 2 x ± d 2 z ± e- lk± - u± 



x 3 2 Xi(y) Ml(u±)] (e^e kl +S^S kl ) 
xdid l z exp{g 4 N cl ( Ul _)x2(y)} ■ 



(55) 



[Aside: At this point, it is easy to verify that the pertur- 
bative result obtained previously in [p|-pyic|] is recovered 
when expanding the exponential to first order. Using 



(e« c W + *yj«) [d x d k zl (u ± )} [did l zl (u x )] 



d 2 p± e 



ip± -u± 



(2n) 2 p\ 



(56) 



the integral over d 2 u± gives (2tt) 2 S(p± + p'j_ — while 
the integral over d 2 b± = d 2 (x± + z±)/2 gives the trans- 
verse area Si. Thus, 



dA 



d 2 k ± dy (2tt) ; 



rtr (|ai 



«2 



2g 6 N c (N 2 - 1) f d 2 p ± 

b± XlX2 



(2 7 r) 3 fc 2 _ 



(2tt) 2 p 2 ± { P± - fcj 



(57) 



This result coincides with those of 0], eq. (36); M, 
eq. (40). The remaining integral has to be regularized by 
introducing a finite color neutralization correlation scale 
A 2 , and can then be written as kj 2 log(fc^/A 2 ) 0. For 



5 



the perturbative regime, that cutoff scale can be chosen 
as A 2 = g^N^/Sn.} 

To simplify eq. ( |55| ) further, note that 

( e tf c « + §V8 kl ) [d* x d k z A] [did[B] = [d 2 x A] [d 2 x B] , 

(58) 

as can be verified most easily in 2-d transverse Fourier 
space: (p±Xq±) 2 +(p_ L -q±) 2 = p 2 L q 2 _(cos 2 (0)+sin 2 (0)) = 
P±1±- F rom the definition of j(u±), see eq. (|45"|), we have 
d 2 ^(uj_) = (2 + log(u 2 L A 2 3CD ))/2vr, and thus 

dN 



/V 2 — 1 



9 Xi(y) 



d 2 /fc ± dy (2 7 r) 4 fc 2 _ 
x (2 + log( M iA| CD )) a 2 exp {g 4 ^c7(«i)X2(2/)} . 

(59) 

We can now integrate by parts. We neglect derivatives 
of the distribution function of the small nucleus, i.e. of 
Xi, and of the logarithm from the propagator. Then, 

^ = V*,(,)|_i/dVdVe-<^ 

x (-2-log(i. 2 ± Aj CD ))exp{ 9 4 Af l7 (n ± ) X2 (!,)} , (60) 

which is just the Fourier transform of 

d^V „ o , ,N 2 - 1 . „ , , o j2 

2ff 2 Xi(2/)^r (- 2 - log(«lA^ CD )) 



d 2 Ui dy * Alvy/ (2tt) 4 
x /d 2 6 ± exp{. 9 4 7V cM ilog( U iA 2 j CD ) X2 (y)/8 7 r} . (61) 



This is our main result. Eq. (|6^) gives the k± and 
y-distribution of produced gluons in the McLerran- 
Venugopalan model, including the renormalization-group 
evolution of x 

In ref. M it was assumed that nucleus 2 represents 
a uniform distribution of charge extending along the 
rapidity axis from yo to y\, such that X2(yo) = 0. 
In other words, neglect QCD evolution of X2 and set 
X^iu) = lAiv ~ Ho)-, with /j,2 = const. Then, integrating 
over rapidity from yo to yi one obtains with logarithmic 
accuracy at x\ <C 1/ Aq CD 



dN N 2 

cx 



d 2 xj_ N c g 2 x 2 j_ 
: j d 2 b ± [1 - exp {-g 4 N c xlt4( yi - y )/8^}] 



(62) 



With /z 2 . = 2ir 2 p re ixG(x) / g 2 (N 2 — 1) one reproduces the 
result of 0]. Here, p re i denotes the (Lorentz-boosted) 
density of nuclcons in nucleus 2: 



Prel 



KR 2 A {yi - y ) 



(63) 



III. DISCUSSION 

In this section we discuss the transverse momentum 
distribution of gluons in various regimes, and the scaling 
of the multiplicity per unit of rapidity with A\ and A 2 . 
When referring to the scaling with the mass numbers 
of the two colliding nuclei, we shall specifically assume 
that at fixed rapidity the color charge densities Xi (y) are 
proportional to A 1 ^ 3 [Q. 

We can understand some general properties of 
eqs. ( |60|j6l| ) even without solving for the RG evo- 
lution of x(y). In the region where £j^qcd <C 
x 2 j_g 4 N c X2(y)/8Tr <C 1, or alternatively Aq CD //c 2 <C 
g 4: N c X2(y)/8Trk 2 _ <C 1, one can expand the exponential 
to first order (the zero'th order term does not contribute 
to k± > 0). Using 



2 - \og(x 2 ± A 2 QCD ) 
j{x±) 



d 2 p ± e ip ^ x ± 
~* pj~ ' 
d 2 q 1 _ e iq± - x ± 
(2^) 2 q\ ' 



(64) 
(65) 



the integral over d 2 u^ in eq. ( |60| ) just gives (2tt) 2 S(pj_ 
q± — k±), and we obtain 



dN 



d 2 6^d 2 A:j_dy 



2g 6 N c (N 2 -l) X i(y)x2(y), k 2 

■ log 



k\ 



g 4 N cX 2/8TT 



(66) 



Thus, one recovers the standard perturbative ~ 1/^1 
behavior at very high fcj_, with a logarithmic correction 
analogous to DGLAP evolution U 10 1. Note that xi, X2 
scale as A 1 ^ 3 and A\^ 3 0, respectively, while the inte- 
gral over d 2 6j_ gives a factor of -kR\ oc A^ 3 . Therefore, 
in this kinematic region dN/d 2 k±dy scales like A 1 / 3 A 2 , 
up to logarithmic corrections. This holds also for the in- 
tegrated distribution dN(k± > po)/dy above some fixed 
^-independent scale po. On the other hand, when inte- 
grating over k\ from g 4: N c X2/8'x to infinity, the contri- 
bution from large k± to the rapidity density is 

dN _g 2 (N 2 -l) 



d 2 6^dy 



-xi(y) ■ 



(67) 



Again, the integral over d 2 6^ gives a factor ttR 2 oc A 2 ^ 3 , 

and so dN/dy scales like A^A 2 / 3 ; in this regard, see 
also the discussion in [jl2j . The transverse energy can be 
obtained from dE± = fc^dA^, using the number distribu- 
tion dfjff 



dE A 



d 2 6^dy 



2AT, 



= ^(N 2 -l)xi(y)Vx2~(yj 



(68) 



Finally, from (|6^) and (^) the average transverse mo- 
mentum in the perturbative regime is 



6 



4Q 



(2) 



(69) 



where Q? ] (y) = Vg 4 N cX 2(y)/8ir. 

Within the "saturation regime" , i.e. when k\ < 
g 4 N c X2{y)/8TT but k\ > g 4 N c xi(y)/&Tr, the upper limit 
on the integral over u± in eq. ( J60| ) is effectively given by 
u\ < 8ir/g 4 N c X2{y)- That is because the exponential 
suppresses contributions from larger itj_. For the deriva- 



tive of the propagator we may again use eq. (f 
we find 



Then 



dA 



d 2 6j_d 2 fc^dy 



2s 2 xi(y) 



A 2 



1 



2wp 2 x 



(2n) 4 

d 2 u ± e- l(k± - p±) - UA 



g 2 xi(y) 



g 2 xi(y) 



A 2 - 1 fdV 
(2tt) 3 



1 



2tt P 2 ± (ki 



N 2 



1 1 



(2tt) 3 k\ 



log 



k 2 



g 4 N cX i(y)/8n ' 
(70) 



This form oc Xl/^± ^ s to be compared with that from 
eq. (|66|), oc XiX2/^jLj valid at high k±. A schematic 
distribution^ in transverse momentum is shown in Fig. ||, 
where Qi stands for \J g 4 N c Xi (y) /Sir. 

Using the DGLAP equation for the transverse evolu- 
tion, we can also express the logarithm times the X x m 
eq. ( |7(j| ) in terms of the unintegrated gluon distribution 
function MM. We write 



xi(y) 



1 



irR\ A 2 - 1 
and move the gluon number 



N?(x,pl) 



(71) 



a s N c 



N?(x,p 2 ± ) 



a — fdx'G{x',pl) 

K Jx 

d 



dlogp^ 



xG(x,p\) 



(72) 



inside the integral over d p± in eq. (|7d). This leads to 



dA 



d 2 b ± d 2 kj_dy irR 2 2vrfc 2 



9 4 JVcXi/8^ 



ttR{ 271-/4 



[xG(x,k 2 ± ) - xG(x,g 4 N c xi/8n)] . (73) 



1 We mention again that we are in fact not able to compute 
the distribution below . It has to be obtained numerically 
using the methods of Q . In Fig. ^| we only express the qual- 
itative expectation that the distribution eventually flattens 
out at very small k± [0,0,0] • 



10" 



H 

Z, 10" 



10 ' 



Q 



d) 



Q 



(2) N 



FIG. 2. Schematic k± distribution for particles produced 
in high-energy p + A collisions (or, more generally, for 
particles produced in A\ + A2 collisions at rapidity y 
such that Xi(j/) X^(v))- l n the perturbative regime, 
dN/dk^dy ~ l/k]_. Inbetween the saturation scales for the 
two sources, dN/dk±dy ~ 1/fex- 

A quantitative computation of the radiation distribu- 
tion requires to determine numerically the CGC density 
scales Xi(y) an d X2(y) from some parametrization of the 
gluon and quark/antiquark distribution functions. Also, 
the fragmentation of the radiated gluons into pions must 
be taken into account. We postpone those issues to a 
future publication. 

From (|70|), the fc^-integrated multiplicity in the 
nonperturbative regime g 4 N cX i(y) / '8tt < k\ < 
g i N cX 2(y)/87T is 



dN 1 2 N 
= 4 9 Xl{y) l^V 2 



d 2 b±dy 



— log 2 ^ 

xny) 



(74) 



Thus, at fixed impact parameter, the multiplicity scales 

1/3 



1/3 

as A x ' , up to the square of a logarithm of (A2/A1) 
For the transverse energy in the saturation regime one 
obtains from (|70|) 



JV« N 2 - 1 



d 2 &j^dy 



2vr (27t) 



xi(y) 



vxi(y) - vxa(y) 



Vx2{y) ]ar X2(y) 
2 xi (2/) 



(75) 



In the saturation regime (0j7|) as well as in the per- 
turbative regime ( |67| , |68| ) the transverse energy per gluon 
is practically independent of A\ 1 while a weak increase 



oc A^ 6 is expected. 



The average transverse momentum in the saturation 
regime follows from ( |74] , |75| ): 



g-l-logg 
log 2 £ 



(76) 



7 



where = yj Xi{y) I X2{y) = Qs \v) I Q K s ''' (y) ■ From 
dimensional considerations, it has been suggested ]l3[j 
that in symmetric A+A collisions, and at central rapidity, 
(k±) 2 scales with the multiplicity per unit of transverse 
area and of rapidity, 



dA 
d 2 6^dy 



(77) 



A similar scaling relation can be derived from eqs. 
for the asymmetric case, 



fa 



dN g 2 (C-l-logO 2 



d 2 b ± dy £ 2 log 6 e 



(78) 



Thus, (k±) 2 is proportional to the multiplicity per unit of 
rapidity and transverse area, times a function of the ratio 
of the saturation momenta. If source one is very much 
weaker than source two, i.e. in the limit | log£| S> 1 — £, 
the third factor on the right-hand-side of ( ffsQ depends 
on log£ only. Neglecting that dependence, and assuming 

1/3 

as before that xi,2 are proportional to A 12 , one has the 
approximate scaling relation 



oc 



Ai 



1/3 



dA 



d 2 b±dy \AiA 2 



1 



1/3 



dN 
dy 



(79) 



In practice though we expect significant corrections to 
the simple scaling relation (|79|), as given by eq. (frsl) . 

Qualitatively, the rapidity distribution predicted from 
eq. ((70|) is as follows^, see Fig. |[ For rapidities far from 
the fragmentation region of the large nucleus, and for 
g^NcXiiy)/^ <k\< g i N cX 2{y)/STT, dN/d 2 k ± dy varies 
like 



d 2 A 
dy 2 



Yi(y) 



(80) 



where we have suppressed the dependence on transverse 
momentum, which is supposed to be held fixed some- 
where within the saturation regime. Thus, an experi- 
mental measure for the RG evolution of the CGC density 
parameter is 



dlogdA/dy dlogxi(y) 



dy 



dy 



(81) 



2 A quantitative computation requires to solve for the RG 
evolution of the x' s first, which is out of the scope of the 
present manuscript. 




FIG. 3. Schematic rapidity distribution for particles pro- 
duced in high-energy p + A collisions (or, more generally, for 
particles produced in A\ + A2 collisions at Ai <C A2). The 
upper curve refers to the perturbative regime, the lower curve 
refers to k± between the saturation scales for the two sources. 

Now consider the case of high fc 2 > g 4 N c X2(y)/8Tr 
described by eq. (|6^). In that regime the rapidity distri- 
bution is proportional to Xi{y)X'2,{y) > and so dA/dyd 2 fci_ 
varies with rapidity like 



dlogdA/dy dlogxi(y) , dlogX2(y) 



dy 



dy 



dy 



(82) 



Subtracting (^TJ) from (g), that is (d 2 N / dy 2 ) / (dN / dy) 
measured at small transverse momentum from that at 
larger transverse momentum, provides an experimental 
measure for the RG evolution of Xi{y)- 



IV. SUMMARY 

In summary, we have computed the radiation field 
produced in the collision of two ultrarelativistic, non- 
abelian, classical color charge sources for the case where 
one of the sources is much stronger than the other. Ac- 
cordingly, we have linearized the Yang-Mills equations 
in field 1, but solved them to all orders in field 2. The 
renormalization-group evolution of the color charge den- 
sity x is not dropped. This problem is relevant for p + A 
collisions at high energy, or more generally for A\ + A2 
nuclear collisions where A\ <C A2, or even for symmetric 
A + A collisions at large rapidities far away from y = 
(i.e. midrapidity) where the RG evolution ensures that 
the effective xi is much smaller than X2- 

We obtain the following results relevant for phe- 
nomenology. At high transverse momenta, k\ > 
g 4 N c X2/8iT, the distribution in k± is proportional to 
the standard XiX2/k 4 L known from perturbation theory. 
Thus, the unintegrated distribution scales like (A1A2) 1 / 3 . 
The total contribution from high transverse momenta, in- 
tegrated over k± and impact parameters b±, scales like 
A\ /3 A 2 2 /3 . 

In the saturation region g 4 A c X2/87r > k\ > 
g 4 N c xi/8TT, the distribution is proportional to Xi/^j.! 
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it decreases much less quickly with transverse momen- 
tum than the result from perturbation theory. This may 
in principle provide experimental information as to the 
value of Xi{v)- At fixed k\, the gluon distribution scales 
like A 1 ^ 3 for fixed impact parameter, or like A\^ 3 A^ 3 
when one integrates over d 2 b±. Note that, up to loga- 
rithmic corrections, the /c^-integrated distribution scales 
in exactly the same way with Ai as in the perturbative 
regime (no matter whether impact parameter selected or 
integrated). In contrast, at fixed k±_ and b± the multi- 

1 /3 

plicity in the perturbative regime scales as A 2 while in 
the saturation regime it is independent of Ai (up to a log- 
arithm) ! (Or, without impact parameter selection, we 
have a scaling with Ai in the perturbative regime versus 
scaling with A^ 3 in the saturation region.) 

Furthermore, at fixed transverse momentum, the quan- 
tity dlog(diV/dy)/dy allows an experimental measure- 
ment of the RG evolution of the color charge density 
parameter x, and a check whether the saturation regime 
has been reached. The slope of dN/dy at rapidities far 
from the fragmentation region of the large nucleus mea- 
sures the RG evolution of xi- Also, subtracting the slope 
of the dN/dy measured at small transverse momentum 
(within the saturation regime) from that at larger trans- 
verse momentum (in the perturbative regime), provides 
experimental access to the RG evolution of X2> and for its 
A dependence. Such differential measurements at RHIC 
and LHC should provide insight regarding high-density 
QCD and the properties of the CGC, for example the 
value of its fundamental parameter x and its RG evolu- 
tion (in rapidity). 
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